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ABSTRACT 


Integral  transforms  are  used  to  solve  singular  initial  value 
problems  for  Stokes'  slow-motion  equations.  The  method  is  applied 
to  investigate  the  decaying  process  of  straight  and  circular  dis¬ 
continuity  lines  as  well  as  the  dissipation  of  local  disturbances 
in  an  infinite  medium.  A  criterion  for  the  occurrence  of  secondary 
vortices  is  derived.  Numerical  results  are  displayed  for  periodic 
initial  disturbances  which  demonstrate  graphically  the  spreading 
of  the  disturbance  from  a  discontinuity  line  into  the  fluid  under 
successive  development  and  decay  of  secondary  vortices.  A  more 
detailed  sequence  of  dissipating  vortices  is  evaluated  numerically 
and  displayed  by  streamline  patterns  in  connection  with  Lamb's 
vortical  eigenmotions  in  an  infinitely  long  cylinder  of  finite 
radius. 
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1.  Introduction 


In  a  recent  paper  [l]  complete  systems  of  time -dependent  sepa¬ 
rable  solutions  of  Stokes1  slow-motion  equations  have  been  constructed 
for  rectangular  and  annular  flow  regions.  This  was  achieved  by  means 
of  a  generalized  separation  technique  introduced  in  [2,  3].  In  addi¬ 
tion,  so-called  associated  separable  solutions  have  been  found  by 
differentiating  or  integrating  separable  solutions  with  respect  to 
an  integration  parameter.  Their  significance  for  the  construction 
of  certain  solutions  of  the  complete  Navier-Stokes  equations  has 
been  pointed  out. 

With  the  aid  of  those  separable  solutions  eigenmotions  have  been 
determined  which  are  required  to  satisfy  certain  incomplete  boundary 
data  without  prescribed  initial  values.  For  given  boundary  and 
initial  data  a  discrete  or  continuous  superposition  of  eigenmotions 
is  necessary.  Because  of  the  complicated  nature  of  this  fitting 
procedure  the  introduction  of  integral  transforms  suggests  itself. 

They  can  be  utilized  because  separable  solutions  exist. 

In  this  paper  the  integral  transform  technique  is  applied  to 
describe  decaying  straight  and  coaxial  discontinuity  lines  as  well 
as  the  obliteration  of  local  disturbances  in  a  laminar  fluid  flow. 

The  analysis  is  restricted  to  plane  motions  so  that  use  can  be  made 
of  the  stream  function  i|r.  Hence,  Stokes1  slow-motion  equations  are 
reduced  to  the  single  fourth  order  partial  differential  equation 

(1) 
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where  A  is  the  Laplacian  operator  and  t  the  product  of  time  and 
constant  kinematic  viscosity.  Numerical  results  are  obtained  for 
selected  examples  which  demonstrate  graphically  the  sequence  of 
decay. 

It  may  be  mentioned  that  under  certain  boundary  and  regularity 
conditions  the  solutions  obtained  can  be  used  to  derive  integrals 
of  the  quasi-linear  Navier-Stokes  equations  in  an  exact  linear 
manner  by  an  iterative  procedure.  This  has  been  shown  in  [l,  2,  3] 
for  various  flow  problems. 

2.  Decay  of  a  Straight  Discontinuity  Line  in  an  Infinite  Region 
A  straight  discontinuity  line  may  be  generated  by  two  uniform 
parallel  flows  of  opposite  direction  at  t  =  0  and  may  be  located 
at  the  axis  y  =  0  in  a  Cartesian  coordinate  system  (x,  y) .  (Fig.  1) . 
In  addition,  a  disturbance  along  the  discontinuity  line  is  super¬ 
posed  which  may  be  of  periodic  or  aperiodic  nature  expressed  by  a 
Fourier  series  or  Fourier  integral,  respectively.  Hence,  a  slow- 
motion  solution  of  Stokes1  equation  (1)  in  Cartesian  coordinates 
i|f  +  2l|f  +  “(♦+♦)*.  (2) 

may  be  specified  by  the  following  initial  and  boundary  conditions. 
The  initial  state  is  described  by  the  basic  flow 

♦b(x,  y>  0)  =  Ae  |y|  ,  (AB  0)  (3) 
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and  the  disturbance 


I|r(x,  y,  o)  = 


+  00 

6(y)  £  A(n)e 

n=  -  oj 
+00 


mx 


^  ®(y)  '  A(n)e^nxdn 


(4) 


The  discrete  or  continuous  function  A(n)  is  assumed  to  yield  con¬ 
vergence  of  the  series  or  integral,  respectively.  Furthermore,- 
A(-n)  =  A(n) ,  where  A  is  the  conjugate  complex  value  of  A.  The 
symbol  6(y)  denotes  Diracfs  delta  function.  The  boundary  conditions 
for  t  ^  0  and  -»  £  x  ^  +  »  are 


y  o  +  oo; 


lim 

|y  H 


1ST. 


-  Ap 


(  +  B)y  “  ±  A0 


(5) 


♦  -  0  , 


(6) 


In  addition  to  the  conditions  (5)  and  (6)  it  is  required  that  i|r 
and  ijr  have  a  sufficiently  monotonic  behavior  for  large  values  of  y. 
For  the  x-direction  it  is  assumed  that  i|/b  anc^  t  are  representable  in 

the  forms 


♦b  (x,  y,  t) 


sB  (y»  t) 


f+  OB 


£  S(y, 


Kx,  y,  t) 


n=-co 


+  C0 


J  s&> 


t,  n)einx 
t,  n)e^nxdn 


(7) 

(8) 
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with  S(-n)  »  S(n).  For  the  sake  of  simplicity  it  may  be  required 
that  the  representations  of  ijr  by  equation  (8)  are  four  times  dif¬ 
ferentiable  in  space  and  once  in  time  under  the  summation  or 
integral  symbol 0  At  this  point  it  may  be  mentioned  that  this 
assumption  restricts  the  manifold  of  solutions  considerably*  For 
instance,  if  general  boundary  data  were  imposed  at  some  lines 
x  =  Xq  and  x  =  xx  ,  it  is  easy  to  show  through  a  finite  Fourier 
transformation  of  equation  (2)  that  the  equation  (11)  below  for 
S(y,  t,  n)  must  be  augmented  by  appropriate  inhomogeneous  terms. 

A  detailed  examination  of  this  deficiency  of  separable  solutions 
will  be  made  in  a  forthcoming  paper  by  Schwiderski  [4]. 

The  initial  and  boundary  data  (3)  through  (8)  form  together 
with  equation  (2)  a  well-set  initial  value  problem.  For  the 
basic  flow  (3),  (5),  and  (7)  the  solution 


with  the  velocity  components 

UB  «  ('l'B)y  =  Ab  erf^-^  ,  VB  =  -  (’lfB)x  =  0  (10) 

is  well  known  in  literature  [5]  and  is  set  aside  for  the  time 
being.  The  remaining  initial  value  problem  for  if  can  be  solved 
in  the  following  way: 
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Insertion  of  the  expansions  (8)  in  equation  (2)  yields  for 


S(y,  t,  n)  the  fourth  order  differential  equation 

bm  -  +  £}yy  +  "a(”3  +  &>  -  9  •  (11) 

By  introducing  the  Laplace  transform 

?(y,  X2-  n8)  =  J*  S(y ,  t,  n)e‘(X2"n2)tdt  (12) 

o 

one  arrives  for  equation  (11)  at 

&  ■  n3)(^  -  9  ■  -  &  -  °* n)  ’  (13) 

where  the  initial  value  S(y,  0,  n)  *  A(n)6(y)#  Equation  (13)  can 
be  integrated  twice  to  yield  the  second  order  differential  equation 


•-  A(n)5(y)  +  cneny  +  dne"ny,  n  /  0 
.  ~  A(o)6(y)  +  c0y  +  d0  ,  n  =  0 


(14) 


The  coefficients  cn  and  dn  are  determined  by  the  boundary  condition 
$(±»,  X8  -  ns)  =  0,  which  one  derives  from  equations  (6)  and  (12). 

It  follows  cn  e  dn  =  0. 

An  attempt  to  Integrate  equation  (14)  for  the  desired  particular 
integral  would  require  a  differentiation  of  equation  (12)  with 
respect  to  y. in  order  to  achieve  the  proper  regularity  at  y  =  0 
for  all  t.  Since  the  delta  function  is  involved  for  t  =  0,  such 
differentiation  of  equation  (12)  cannot  be  carried  out  under  the 
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integral  sign.  This  difficulty  is  avoided  by  applying  to  equation 

(14)  the  Fourier  transform 

+00 


fn,\(p)  =  J*  -  n2)eipydy  .  (15) 

9  -00 


Then,  one  obtains 

f  -  — 4.Cn) 

n,\(p)  x3  +  p2 


(16) 


In  [6,  7]  the  original  function  of  fn  ^(p)  with  respect  to  y  is 
tabulated  and  becomes 


$  .  A(n)  e“X  !y  I 

2\ 


(17) 


This  again  is  the  image  function  of  S  and  can  be  inverted  to 


s  =  7^-)exp  •(nS|:  +  fe)  • 


Zv/nt 


(18) 


Integral  (18)  belongs  to  the  class  of  associated  separable  solu¬ 
tions  which  have  been  found  in  [1]  .  Substitution  of  (18)  in  (4) 
renders  the  solution  of  the  initial  value  problem  (2) ,  (4) ,  (6) , 
and  (8) : 

inx  -  n2fc 

(19a) 


y2 


♦<*,  y,  t)  =  — ^e" 
VTTt 


"V°  A  t  \  inX  “  n2  t 
E  A(n)e 

n=-oa 


/“  A(n)einx  “  n@tdn 


(19b) 


Solution  (19a)  describes  the  decay  of  a  periodic  disturbance 
which  is  confined  at  t  =  0  to  the  single  line  y  =  0  and  which  spreads 
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with  advancing  time  into  the  fluid  in  form  of  a  discrete  spectrum 
of  vortical  motions.  Each  mode  n  consists  of  2n  vortices  in  the  interval 
from  0  to  2tt.  For  inodes  n  >  0  the  dissipation  of  the  motions  (19a)  is 
essentially  governed  by  the  exponential  function  exp(-n2t).  Since  its 
damping  coefficient  is  n3 ,  vortices  of  high  order  n  decay  faster  than 
those  of  low  order  n.  Furthermore,  those  vortices  with  largest  coef¬ 
ficient  A(n)  dominate  at  the  beginning  and  become  visible  as  "secondary 
vortices.11  The  2n  vortices  of  this  dominating  mode  grow  in  their 
spatial  extent  and  die  away  as  soon  as  they  are  overcome  by  vortices 
of  lower  modes.  This  mechanism  explains  the  observed  occurrence 
of  small  vortices  of  large  number  in  the  beginning  of  the  decaying 
process  of  a  discontinuity  line  and  their  replacement  by  fewer  vortices 
at  a  later  time  (see  the  experiments  by  Weske  and  Rankin  [8]). 

The  process  described  above  may  be  illustrated  by  a  series  of 

flow  patterns  at  various  times.  In  Figs.  2  through  13  computer 

output  is  displayed  for  streamlines  which  differ  from  each  other 

by  an  equal  amount.  For  convenience  the  Fourier  coefficients  A(n) 

may  be  expressed  by  the  amplitude  Cn  and  the  phase  angle  such 
icvn 

that  A(n)  =  Cne  .  In  the  first  sequence  of  figures  the  motion 
consists  of  the  superposed  modes  n  =  1  and  n  -  2.  The  ratio  of 
their  amplitudes  is  C2  :  Cx  =  10  so  that  the  four  vortices  of  n  =  2 
dominate  at  the  beginning  (Fig.  2).  Furthermore,  the  phase  angles 
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of  the  two  modes  may  differ  by  tt/4  .  If  time  passes  by,  the  second 
vortex  diminishes  and  its  two  neighboring  vortices  join  each  other, 
whereas  the  fourth  vortex  gains  considerable  strength  (Figs.  3 
through  6).  The  remaining  two  vortices,  then,  approach  the  state 
which  is  represented  by  n  =  1  (Fig.  7)  and  finally  die  away  when 
t  tends  to  infinity.  The  same  vortex  configuration  but  with  equal 
phase  angles  (Ca  :  Cj  =  10 ,  as  -  =  0)  dissolves  in  a  symmetric 

way  by  supporting  and  suppressing  two  adjacent  vortices  each  (Figs. 
8  and  9).  In  a  third  sequence  of  figures  a  disturbance  of  mode 
n  =  2  is  superposed  on  the  basic  flow  (9)  with  Cs  :  A0  =  5a/tT"  and 
aa  =  tt/2  .  Fig.  10  shows  the  flow  field  in  an  early  stage.  With 
progressing  time  those  vortices  are  weakened  which  have  flow  direc¬ 
tions  opposite  to  the  main  stream  (Fig.  11).  They  finally  vanish 
(Fig.  12),  and  the  remaining  two  vortices  are  gradually  absorbed  by 
the  basic  motion  (Fig.  13). 

From  solution  (19a)  one  can  derive  the  criterion  for  the 
generation  of  any  secondary  vortex  which  accompanies  a  dissipating 
discontinuity  line:  There  must  exist  at  least  one  mode  n  ^  0  with 
nonvanishing  A(n)  which  is  sufficiently  strong  in  comparison  with 
A(o)  and  AB . 

It  may  be  noted  that  the  solution  (19a)  exhibits  the  same  decay 
ing  mechanism  as  a  complete  solution  of  the  Navier-Stokes  equations 
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As  was  shown  in  [l]  for  the  similar  case  of  Lambfs  eigenmotions 
(see  section  4),  the  solutions  of  Stokes1  and  Navier-Stokes *  equa¬ 
tions  approach  each  other  during  the  final  phase  of  decay  and 
become  equal  if  n  -  0  prevails. 

The  decay  of  an  aperiodic  disturbance  may  be  illustrated  by 
the  example 


A(n)  ■  1,  J 


einxdn 


=  2tt  6  (x) 


(20) 


Then,  the  evaluation  of  the  Fourier  integral  in  (19b)  yields  the 
solution 


i(f  =  exp  -  ^(x3  +  ya)  . 


(21) 


Equation  (21)  describes  the  decay  of  a  single  vortex,  where  the 
stream  function  depends  only  on  the  radial  distance  r  *  (x2  +  y2)^ 
from  the  center  of  the  disturbance.  This  leads  to  the  study  of 
disturbances  in  a  polar  coordinate  system,  which  will  be  the 
subject  of  the  next  section. 

3.  Decay  of  a  Coaxial  Discontinuity  Line  in  an  Infinite  Region 

In  literature  the  solution  of  a  decaying  coaxial  discontinuity 
line  seems  not  to  be  known  which  is  produced  by  a  jump  in  the 
angular  velocity  of  a  fluid  rotating  like  a  solid  body.  This 
problem  corresponds  to  that  of  equation  (3)  for  two  parallel  flows 
of  opposite  direction  and  is  not  considered  here.  Rather  the 
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investigation  is  focused  on  disturbances  as  in  the  previous  section. 
Different  from  equation  (4)  univalent  disturbances  of  a  coaxial  dis¬ 
continuity  line  can  only  be  of  periodic  nature.  Therefore,  if  r  =*  a 
designates  the  location  of  the  singular  line,  the  disturbance  is 
assumed  to  be  of  the  general  form 

00 

t(r,  cp,  0)  «  6(r  -  a)  S  ^cosCncp^-  an)  9  (22) 

n=o 

For  all  t  ^  0,  the  disturbance  shall  vanish  at  r  ■  «,  that  is 
r  «=  a  :  f  =  0  ,  =  0  0  (23) 

The  solution  can  be  sought  as 
00 

Kr,  cp,  t)  =  E  G„(r,  t)cos(ncp  +  an )  ,  (24) 

n=o 

which  reduces  Stokes'  slow-motion  equation  (1)  in  polar  coordinates 
(r,  cp)  to 


(f  +1  a  -  n2  W  +  1  b_.IL. 
Var8  r  3r  r^Aar^  r  Sr  r2 


ia_ 

r  3r 


,0) 

(27) 
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where  G^Cr,  0)  =  An§(r  -  a).  Double  integration  of  (27)  yields 
the  second  order  differential  equation 


1  A_ 

r  dr 


-An6(r-a)  +  c^11  +  dnr"n,  n  ^  0 


-A06(r-a)  +  c0log  r  -t-  d0,  n  =  0 


(28) 


The  boundary  data  (23)  and  the  regularity  condition  at  r  =  0 
require  =  dn  =  0  for  a  >  0.  In  order  to  integrate  equation  (28) 
the  Hankel  transform  of  order  n  must  be  applied, 

00 

hn  X(p)  =  J  ^  5n(r)Jn(rpVrp  dr,  p  >  0  ,  (29) 

*  o 

where  Jn  is  the  n  order  Bessel  function.  This  transformation 
leads  to 

00 

(p3  +  XS)hn  ^  “  An  J*  6(r  "  a)Jn(rp)Vrp  dr  (30) 

or  by  evaluating  the  integral  pn  the  right  to 

hn,7i  =  ps  +  p-DWp  Jn<aP)^  >  a  >  0  ,  (31) 

This  function  is  readily  inverted  with  respect  to  t  and  yields 

H  (p)  =  A  aVp  J  (ap)e"p  fc  .  (32) 

n  n  n 

According  to  [6,  Vol.  II,  p.  51]  the  original  function  of  H^Cp)  is 

Gn'r>  «  "  £  V*n(g'  >*>(-  <33) 
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with  In  as  n™  order  modified  Bessel  function  of  first  kind.  Hence, 


i|f(r,y,t)  =  ^exp(“  AnIn(|^)cos(ncP  +  ^  >  a>0  •  C34) 

-£ 

Since  the  asymptotic  behavior  of  In^ft^^ft)  exp(ft)  *'s  :‘-nc^ePen^“ 
ent  of  n,  the  initial  state  (22)  and  the  boundary  conditions  (23) 
can  easily  be  verified.  The  motion  which  is  represented  by  the 
integral  (34)  corresponds  to  that  of  (19a)  for  a  straight  discontinuity 
line.  Therefore,  its  description  is  here  omitted. 

It  may  be  mentioned  that  the  special  case  n  ~  0  of  the  integral 
(34)  appears  in  [9]  as  the  vorticity  ou  of  a  decaying  discontinuity 
line.  ^  Since  this  solution  fulfills  the  initial  condition 
cu(r,  cp,  0)  «  6(r  -  a),  the  discontinuity  line  considered  in  [9]  is  of 


the  form 


u(r,  'Cp,  0)  s  0,  v(r,  cp,  0)  = 


—  for  r  ^  a 
r 

0  for  r  <  a 


where  u)r  =  (rv)_  • 
r 

A  slightly  different  situation  arises  for  a  =  0  if  the  coaxial 
discontinuity  line  degenerates  to  a  local  disturbance.  Then,  the 
integral  to  the  right  of  equation  (30)  is  zero.  However,  if  one 
changes  the  coefficients  in  (22)  from  to  A^r"1^11  one  arrives  at 


the  image  function 

,  p*+£ 

h  .  =  A  2’n  a  TTs 
n,A  n:  n  p  +  A 
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whereby  the  coefficient  dn  in  equation  (28)  must  vanish  in  order 
to  obtain  a  finite  expression  for  According  to  [6,  Vol.'ll, 

p.  63]  the  inversion  of  (35)  with  respect  to  r  results  in 


*n  -  £  • 

i.  t  - 

where  is  the  ntn  order  modified  Bessei  function  of  second  kind 
Finally,  the  inverted  function  of  (36)  with  respect  to  t  is 


Gn(r,t) 


.n 


(2t) 


n+i 


5 


(36) 


(37) 


(see  [7,  p.  128]).  Formally,  the  integral  (37)  can  be  obtained  from 
(33)  if  one  replaces  A^  by  Ana"1P"n  and  determines  the  limit  for 
a  ->  0. 

The  integral  (37)  belongs  to  the  class  of  associated  separable 
solutions  which  have  been  constructed  and  discussed  in  [l].  They 
are  interpreted  as  "explosive  rotations.11  The  special  case  n  =  0 
may  be  compared  with  equation  (21).  It  follows  that  the  initial 
conditions  i|/(x,  y,  0)  =  2TT6(x)6(y)  and  \|/(r,  y,  0)  *»  ~  6(r)  are 
equivalent.  Indeed,  they  lead  to  the  same  solution. 

4.  Decay  of  Vortices  in  an  Infinitely  Long  Cylinder  of  Finite  Radius 

The  sequence  of  decaying  vortices  from  higher  modes  downward  to 
lower  ones,  which  have  been  described  in  section  2  for  a  straight 
discontinuity  line  in  an  infinite  flow  region,  shall  now  be  demon- 
strated  by  means  of  a  vortex  configuration  in  a  finite  circle, 
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since  this  process  has  much  resemblance  to  the  phenomenon  observed 
by  Weske  and  Rankin  [8].  The  vortices  at  t  =  0  may  be  considered 
as  an  intermediate  state  of  a  dissipating  discontinuity  line.  As 
was  pointed  out  in  [l]  such  a  decay  can  be  described  by  Lamb’s 
eigenfunctions 

'Wr  W08*”*  +  Qmn>  tJn(Xmnr)-Jn(\nn)rn>xPC-4n  t}  08) 

Jn+i  (^mn)  “0  n=0,  1,  2.*.,  m=l52,3,  ...  (39) 

where  the  second  equation  determines  the  eigenvalues  and  where 

the  radius  of  the  cylinder  is  unity. 

For  numerical  calculations  a  vortex  ensemble  has  been  selected 
that  consists  of  the  three  modes  m  =  1;  n  =  2>  3,  4.  They  represent 
4,  6,  and  8  vortices  each.  The  amplitudes  are  arranged  to  have  the 
ratios  A4:A3:AS  =  10,000:1000:1  so  that  the  mode  n  -  4  dominates  at 
t  =  0  ov£r  the  other  two  (Fig.  14),  and  so  that  n  =  3  prevails  over 
n  °  2  at  a  later  time.  The  phase  angles  are  zero.  In  order  to 
exhibit  the  essential  features  in  the  decaying  process  the  follow¬ 
ing  steps  are  selected:  At  t  =  0.11  (Fig.  15)  the  flow  field  in 
the  lower  half  of  the  cylinder  has  weakened  in  such  a  way  that  the 
vortex  which  was  near  the  180°  mark  in  Fig.  14,  does  not  appear  any 
more,  and  that  the  two  neighboring  vortices  are  about  to  join.  The 
latter  incidence  happens  at  t  =  0.13  (Fig.  16).  Six  vortices  are 
now  visible  and  regroup  until  an  almost  symmetric  configuration 
(that  is  the  prevailing  mode  n  =  3)  is  reached  at  t  =  0.25  (Fig.  17). 
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Then,  the  same  process  repeats  itself  and  ends  with  four  dominating 
vortices  (Figs.  18  through  21)  which  dissipate  slowly  when  t  -*  ®  . 
The  asymmetric  decay  is  owing  to  the  property  of  the  trigonometric 
function  and  can  be  observed  in  the  photographs  by  Weske  and  Rankin 
[8]. 

It  may  be  mentioned  that  the  period  of  time  from  one  dominating 
mode  to  the  other  can  be  minimized  by  selecting  optimal  values 
for  the  amplitudes  and  phase  angles.  This  time  interval  increases 
with  lower  modes.  For  the  example  above  the  ratio  of  the  time 
intervals,  which  are  required  to  expose  the  vortices  of  modes  n  =  3 
and  2  respectively,  is  1:2.2.  This  compares  favorably  with  Weske 
and  Rankin's  observation  of  1:2.5 '[10 1  . 
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Fig.  3:  Streamlines  of  a  decaying  straight  discontinuity  line  with  C; 
01  p  “  Q?t  =  tt/4  at  t  “  0  a  6  o 


Fig.  5:  Streamlines  of  a  decaying  straight  discontinuity  line  with 
Ca  :  Ci  =  10,  aa  -  or,  =  tt/4  at  t  =  1.6  . 


Streamlines  of  a  decaying  straight  discontinuity  line 
with  Cs  :  Cn  =  10,  a2  -  os  =  tt/4  at  t  =  2#0  . 
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Fig*  11:  Streamlines  of  a  decaying  straight  discontinuity  line 
with  C2  :  Ab  =  5\At  ,  a2  =  tt/2  at  t  =  0,2 
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